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By using a propagation scheme for current matrices and an auxiliary mode expansion method, we 
investigate the transient dynamics of a single molecular junction coupled with a vibrational mode. 

Our approach is based on the Anderson-Holstein model and the dressed tunneling approximation 
for the electronic self-energy in the polaronic regime. The time-dependent currents after a sudden 
switching on the tunneling to leads and an abrupt upward step bias pulse are calculated. We show 
that the strong electron-phonon interaction greatly influences the nonlinear response properties of 
the system, and gives rise to interesting characteristics on the time traces of transient currents. 
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I. INTRODUCTION 

The charge transport through nanoscale physical sys¬ 
tems, including semiconductor quantum dots,^ single 
molecular junctions(SMJ),^ and carbon nanotubes,etc., 
is currently under active investigations in the field of na¬ 
noelectronics. During the electron tunneling processes, 
the charging of a molecule or a nanotube quantum dot 
may cause deformations of its surrounding geometry and 
the emission or absorption of phonon(vibron) quanta, 
inelastic electron tunneling effects due to the electron- 
phonon interaction (EPI) are ubiquitously observed.^”® 
It was found that the EPI can give rise to a variety of 
physically interesting phenomena in experiments on the 
dc characteristics of SMJs or the suspended carbon nan¬ 
otube quantum dots, e.g., strong current suppression at 
low bias voltage termed Eranck-Condon blockade'^, side 
peaks of differential conductance®, and negative differen¬ 
tial conductances in the current-voltage characteristics®. 

The dynamical properties of nanoscale systems under 
time-dependent external potentials are also of great in¬ 
terest and importance for their potential application in 
electronics devices. Some time-dependent potential re¬ 
lated phenomena have been investigated for the systems 
in the linear response and the nonlinear regimes, for in¬ 
stance, the ac response and current noise spectra,the 
quantum charge pumping^^“^® and the photon-assisted 
tunneling^d-i^ The theoretical investigation of time- 
dependent transport or transient properties of nanoscale 
electronic systems is a challenging task, in particular for 
the systems with strong EPI or strong Coulomb interac¬ 
tion. For a noninteracting quantum dot, a general the¬ 
oretical framework for time-dependent transport proper¬ 
ties has been given in the seminar work by Jauho et 
and the switching effect of currents in response to sharp 
step- and also square-shaped voltage pulses has been an¬ 
alytically obtained by Maciejko et al}^ . The transient 
current through a quantum dot with weak Coulomb in¬ 
teraction were studied by Schmidt et al}'^ within the 
Anderson impurity model by using a mean field the¬ 


ory approximation and also the nonequilibrium Monte 
Carlo simulation method. The transient and ac response 
properties of the Anderson impurity model in the Kondo 
regime were addressed by using the time-dependent non¬ 
crossing approximation. Recently, we have studied 
the time-dependent electron transport through a quan¬ 
tum dot with strong on-site Coulomb interaction^ * using 
a second-order quantum rate equation and a propagation 
scheme proposed by Croy and Saalmann.^^ 

For the system with EPI, e.g., a SMJ or a carbon 
nanotube quantum dot coupled with vibrational mode, 
an archetype theoretical model is the Anderson-Holstein 
model. The transient effects of the Anderson-Holstein 
model in the weak EPI limit have been investigated using 
perturbation theory.^® For a fast phonon mode (wq ^ P) 
and in the case of sudden switching on the coupling to the 
leads, it was shown that the time-dependent total current 
exhibits plateau structures, which is attributed to the ef¬ 
fect of electron shuttling due to the EPI. The transient 
dynamics of the Anderson-Holstein model in the strong 
EPI region is more complicated, a variety of theoretical 
methods were used to address this problem, e.g., nonequi¬ 
librium Green’s function technique,the polaron tun¬ 
neling approximation,^® the diagrammatic Monte Carlo 
simulation approach,^®”®® multilayer multiconconfigura- 
tion time-dependent Hartree method®®”®®, etc. Bistabil¬ 
ity of steady state currents and the dot occupation num¬ 
ber in the strong EPI case were shown for different initial 
preparations of the electronic or the phonon states of the 
system.^®’®®’®®”®® However, a general consensus about the 
dynamics in the strong coupling region and the condition 
for bistability in the Anderson-Holstein model haven’t 
been obtained.®® 

In this paper we will study the transient dynamics 
of the Anderson-Holstein model in the strong EPI re¬ 
gion by using the Lang-Firsov unitary transformation 
and a recently proposed approximation for the electron 
self-energy termed the dressed tunneling approximation 
(DTA).®^’®® It was shown®® that the DTA eliminates 
some pathologies of the single-particle approximation 
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(SPA) or the polaron tunneling approximation (PTA), 
and provides results for the current and the differential 
conductance in good agreement with that of the more 
elaborated methods such as the diagrammatic Monte 
Carlo simulation.^® We consider the system in the pola- 
ronic regime with a fast phonon mode, and show that the 
time-dependent total current and also the displacement 
current exhibit obvious signatures of EPI effect when the 
system is under a sudden switching on of the tunnel cou¬ 
pling to leads or an upward step-shaped bias voltage. 
Therefore, the measuring of transient currents can give 
us important information on the vibrational modes in the 
SMJs. 

This paper is organized as follows: In Sec. II the model 
Hamiltonian and the self-energy terms in the presence of 
external time-dependent potentials are given. In Sec. Ill 
within the propagation scheme and the auxiliary-mode 
expansion method, we derive equations of motion for the 
dot occupation and the auxiliary current matrices. In 
Sec. IV numerical results and discussions on the transient 
currents are presented. In Sec. V the summary and 
conclusion are given. 


II. MODEL HAMILTONIAN AND THE 
SELF-ENERGY TERMS 


The Anderson-Holstein model can describe the elec¬ 
tron transport between a molecular quantum dot(QD) 
and two metallic leads, where it is assumed that the 
molecular QD has only one energy level (with energy ed) 
involved in the electron tunneling process, and is also 
coupled to a single vibrational mode (phonon) of the 
molecular with the frequency ujq and the EPI strength 
gep. The Hamiltonian of this model is given as 


H ^ ekrj{t)c\^Ckp + ed{t)d''d + ujoa^a + gepd^d{a) -I- a) 

krj 


krj 


( 1 ) 


where r] = L,R denotes the left and right leads, €kp{t) = 
ekr, + \it) and edit) = -b Aq( t), with A,,(t) and Aq( t) 
being the time-dependent potentials applied to the lead 77 
and the dot region, respectively, jp (t) describes the time- 
dependent tunneling-coupling matrix element between 
the QD and lead 77, and takes the form of 'jpit) = 'ypu{t), 
with uit) being a time dependent function, e.g., in the 
case of sudden switching on the tunneling at t = 0, 
uit) = 0 (t), with 0 (t) denoting the Heaviside step func¬ 
tion. We can also define coupling strength between the 
lead and the QD as r7,(e) = 27 r^^ \ jp\^Sie — e^). It 
should be noted that here we neglect the on-site Coulomb 
interaction in the molecular QD and consider a spin¬ 
less electron model. Applying the Lang-Firsov unitary 
transformation:®® H = e^He~^, with S = gd^dia^ — a) 
and the dimensionless parameter g = gepjujQ. The trans¬ 


formed Hamiltonian becomes 

H = E ^kr,it)c\^Ckp -b edit)dU + ujoa'^a 

krj 

+ E \lvit)Ar,dX + Ckp 

kr} 


( 2 ) 


where the phonon shift operator X = = 

g-g(a-at)^ and edit) = ed + Ao(t) with the renormal¬ 
ized energy level ed = ed — g'^pjujo. In the transformed 
Hamiltonian, the direct coupling between the electron 
and the phonon is eliminated, but the dot-lead tunneling 
amplitude is modified by the phonon shift operator A, 
which is responsible for the observation of the Franck- 
Condon steps in the dc current-voltage characteristics of 
the SMJs.^ 

The electric current from the lead 77 to the molecular 
dot Ipit) = -e(^), is given by 


^ ^E [^vit){(iridX)--f*it){XU'^Ckp) 

k 

= |[n,(t) + n;(t)], ( 3 ) 

where a current matrix n^(t) = i J^k lriit){c\pdX) is in¬ 
troduced. By using the equation of motion method, one 
can express the current matrix as an integration on the 
closed-time contour over the combination of the GFs of 
QD operator and its corresponding self-energy as follows 

Upit) = f dt'[G^(t,t')S<(t',t) + G<(t,t')S“(t',t)] 

J —oo 

= f dt'[G>it,t')i:<it',t)-G<it,t')j:>it',t)]. 

J —oo 

( 4 ) 

Here Gc(t,t') = (Tcd(t)d^(t')) denotes the GF of dot 
operator in the transformed Hamiltonian of Eq. ( 2 ), and 
satisfies the Dyson equation on the closed-time contour 


[i^^-edit)]Gcit,t') = 5it,t')+^ J dtiT.pit,ti)Gciti,t') ■ 

( 5 ) 

The self-energy term E^(t',t) is given by 


= '^-f*it’)gkpit',t)-fpit)KitG') , (6) 

k 


where gkpit',t) is the GF of lead 77 in the presence 
of the time-dependent potential A,,(t), and A(t,t') is 
the propagator of the phonon shift operator, A(t,t') = 
(TcA(t)Al(t')). By making comparison with the nonin¬ 
teracting resonant tunneling model, one can see that here 
the self-energy due to tunneling into the leads is dressed 
by a polaronic cloud term. 

In this work we consider the system with a fast phonon 
mode (wq ^ T), and a strong energy dissipation of the 
vibration mode to a thermal bath, e.g., a substrate or a 
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backgate. Then we can assume the phonons are always 
in equilibrium and have an equilibrium Bose distribution 
ub = — 1) at the temperature T = 1//3. The 

phonon shift operator GF K{t, t') can be obtained by its 
equilibrium correlation function,and its greater and 
lesser parts are given explicitly as 

= exp{-g^[nB{l - 

+ + (7) 

and They can also expressed as 

sums of Fourier series 


n— — oo 
oo 

K<{t,t')= . (8) 

n— — oo 

The factor Wn is the weighting factor describing the elec¬ 
tronic tunneling involving absorption or emission of n 
phonons. 

Wn = e-9^(2"B+i)j^(2gVnB(^B + , (9) 

where In{x) is the nth Bessel function of complex argu¬ 
ment. Substitute Eq. ( 8 ) into the self-energy expression 
Eq. ( 6 ), we obtain the greater and lesser parts of self¬ 
energy as follows 

k 

oo 

= —iu(t')u(t) Wn 

n— — oc) 

drA,(r) 


de 

27T 


[l-/.(e)] 


k 

oo 

= iu{t')u{t) W-n 

n— — oc 

drA,(r) 

Here / 17 (e) = -I-1] is the Fermi distribution 

function, with being the chemical potential of lead 77 . 



III. AUXILIARY-MODE EXPANSION AND THE 
PROPAGATION SCHEME 

In the above section we show that an energy integral 
over Fermi function is involved in the calculation of the 
lesser or greater self-energy. It will be a heavy burden 
on computation if multiple integral is encountered in the 
study of transient dynamics. At finite temperature, a 
well-known method for performing the energy integral 


is to use the reside theorem and represent the integral 
by the Matsubara summation,"^'* but this method is suf¬ 
fered from slow convergence especially at low tempera¬ 
ture. Some highly accurate and more efficient methods 
have been developed in the literatures. In this work, 
we will use the auxiliary-mode expansion introduced by 
Croy and Saalman.^^ The Fermi function is expanded in 
a sum over simple poles 


/^(e) 


1 lyA/ I ^ 1 \ 

2 ve-xi e-Xw/ 


( 11 ) 


where x^p = Mi? i *Xp//S (with Xp > 0)- Thus all poles 
Xpp iXpp) s-re in the upper (lower) complex plane (details 
for how to calculate these pole parameters can be found 
in the appendix of Ref. [22]). Np is the number of pole 
pairs. Therefore, we will replace the Fermi function in 
Eq. (10) by the auxiliary-mode expansion above. 

Next, in order to take into account finite-bandwidth 
effects, we take the energy dependence of the linewidth 
function r^(e) to be Lorentzian 


^,{e) 


e^ + W^ 


r+ 

1 ? 

e-iW 


iW 


( 12 ) 


where r° is the linewidth constant and W is the band¬ 
width. It indicates that this linewidth function has two 
simple poles at e = iLiW with residues F^J^ = /2. 

Then Eqs. (11) and (12) are plugged into the definition 
of the self-energies in Eq. (10), and the energy integral is 
evaluated by the contour integration method. Eor t > 
we close the contour in the upper half plane and pick up 
the poles to obtain 


S>(t',t) = u{t')u{t)Ywn\ r+[l - f^ilW)] 


Nf ^ 

P=1 1 


dTA,(T) 


(13a) 


Y = -u{t')u{t) Y w-n jr+ {iW) 

n ^ 

Nf 

-t) _ -t) 

p^l 

.g-*//drA,(r) ^ (13b) 

where {iW) denotes that the Eermi function at the 
pole iW calculated by the expansion given in Eq. (11). 
To simplify the notation, one can introduce the coeffi¬ 
cients and exponents as follows 

= { ± WnT^il - {±zW)],±^WnTr,ixt)] > 

(14a) 








4 


= f^TW-nr^f^{±lWl±^W-nT,ix^p)^ , (14b) 
X^nX = \ ±iW + nwo, X% + nwo L (14c) 


The Eqs. (20), (21) and (22) form a closed set of first 
order differential equations for current matrices in this 
auxiliary-mode expansion method, and give the basis of 
a propagation scheme for the transient dynamics of this 
EPI system. 


Then the self-energies can be written in a compact form 
E|(t',f)=,r(t)^Ej„^(f',t), (15) 


nX 


in which for t > t', 

sJ„A(i'>0 , (16) 

and for t < t', we have 


The current matrix 11,, (t) in Eq. (4) can be written 
in terms of auxiliary-mode expansion current matrices 
n?,nA (l) 

n,,(f) = u{t) ^ n,,„A(i), (18) 

nX 


with 


n,nA(t)= f dt'[G>(t,t')S<„;,(t',t)-G<(t,f')S>„;,(t',t)] 

J —CO 

(19) 

By differential the above equation with respect to the 
time variable t, and follow the same procedure as in Ref. 
[22], we obtain the equation of motion for the auxiliary 
current matrices 


= u{t) T"<’+ 


r^nt + ^d(0(d-r<;+) 
ed-f Ao(t)-(x+„^-l-A,,(f)) n,,„A(t) 


+ u{t) ^ ^r]nX,r]'n'X': (^9) 

ri'n' X' 


where a new quantity ^T^nX.rj’n’X’it) is introduced, which 
satisfies the following equation of motion 




^rjnXjTf'n'X' (i) — ^(^) 




Tj'n'X' 


+ (Knt-K^xK'n'X'(t) 


(xU + \(t)) 


(Xn'n'X' + 


^r]nX,r]'n'X' (0 • 


( 21 ) 


Erom the current conservation condition in the QD re¬ 
gion, it is easy to see that equation of motion for the QD 
occupation number nd{t) is given by 

|nd(t) = 5^[n,(f) + n;(t)] 

7 ] 

= u{t)J2[^vnx{t) + n;^^{t)]. (22) 

rjnX 



FIG. 1. Left, right and total current lL{t), 7fl(t) and I{t) 
as a function of time at different electron-phonon coupling 
constants: g = 1.0 (a); 2.0 (b), respectively. The energy 
parameters used are in units of F, with the bandwidth of 
leads W = 50, F® = F^ = F = 1.0, the bias voltage A/r = 2.0, 
the frequency of vibration mode wo = 10.0, and the dot level 
Ed = —5.0. The quantum dot is assumed initially in the empty 
state. 


IV. RESULTS AND DISCUSSIONS 

By solving the above set of partial differential equa¬ 
tions with the finite difference method, we can obtain 
the numerical results of the transient current for vari¬ 
ous initial conditions. For simplicity, we consider the 
system with symmetric tunnelling coupling to the leads, 
with r° = = 1.0, and T = (T^ + r5j)/2 is taken as 

the unit of energy. The phonon frequency is assumed as 
Wo = 10.0, and the Fermi level of the leads at equilib¬ 
rium give the reference point of energy, (xl = TI'R = 0. 
In our numerical calculation, the system is taken as at 
a finite temperature T = 0.5. The number of poles 
Np = 80 is used for the expansion of Fermi distribution 
function, which ensures the convergence of the numerical 
results. We study the transient properties of the cur¬ 
rent in this system, which can be divided into two terms: 
the total current I{t) = [lL{t) — /fi(t)]/2, which mea- 
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sures the charge transported through the system, and 
the displacement current I dispit) = I hit) + Init), which 
reflects the change of occupation number in the dot re¬ 
gion Idispit) = nd{t)p-'^ 



FIG. 2. Time-dependent total current I{t) for initially empty 
molecular quantum dot. The electron-phonon coupling con¬ 
stants g = 2.0, the bias voltage Ap = 2.0. The other param¬ 
eters used are the same as that in figure 1. 



t(1/r) 


FIG. 3. The displacement current Idispit) for initially empty 
molecular quantum dot. The electron-phonon coupling con¬ 
stants g — 2.0, the bias voltage A/i = 2.0. 


A. Sudden switching on of the coupling to leads 

We first consider transient current behaviors when the 
coupling of the quantum dot to leads is suddenly switched 
on at an initial time {to =0). In this case the time- 
dependent function u{t) = 9{t), and we will assume the 
quantum dot initially to be empty. In Fig. 1 the cur¬ 
rent II and Ir from the left and the right leads and also 
the total current I{t) are plotted for the system with a 
bias voltage A/i = 2.0. Because the finite band-width 
(we take W = 50) of leads has been taken into account 
in our calculation, the current Ir and Ir are zero at the 
initial time t = 0, then increase to the maximal values at 
a very short time-scale of f ~ XjW. It is well known that 
in the infinitely wide band approximation, an unphys¬ 
ical current jump to a finite value immediately at the 
initial time is obtained. Therefore, the finite bandwidth 
approach adopted in this work gives a more accurate de¬ 
scription of the transient current behavior. The time 
evolution of currents in Fig. I shows that the system 
will reach to a steady state in the long time limit, how¬ 
ever a much longer time is needed for the system with the 
strong FPL One can attribute this phenomenon to the 
strong renormalization of the tunneling coupling between 
the leads and the quantum dot because of the EPI. For 
weak EPI case {g = 1.0) shown in Fig. 1(a), the time- 
dependence of currents is relatively smooth. However, 
for the strong EPI case {g = 2.0) in Fig. 1(b), more com¬ 
plex structures of current evolution are clearly indicated 
in the short time region. 

In Fig. 2 we examine the total current in the short 
time region for the systems with strong EPI in more de¬ 
tails. Three different dot level positions are considered 
for the SMJ under a fixed bias voltage (A/r = 2.0). As 
the dot level = 0, which is located in the transport 
window {ijlr < < /^i), the steady total current in the 

long time limit has a relatively high value. On the con¬ 
trary, the steady total current is very small when the 
dot level td is far below or above the Fermi levels of the 
leads. It shows that in the short time region the total 
current exhibits Fano lineshape in the time domain with 
peak and dip structures at the time t k. ‘I'KnjijJo-, where n 
is an integer. Similar current oscillating behaviors were 
also found previously by the time-dependent APTA and 
the diagrammatic Monte-Carlo simulation. and were 
attributed to the electron shuttling due to the EPI. This 
result indicates that experimental measuring of the tran¬ 
sient total current through SMJs can provide useful in¬ 
formation about the vibration modes of the molecule be¬ 
tween the leads. 

The displacement current Idispit) which reflects the 
charge accumulation process in the dot is plotted in Fig. 
3. The displacement currents for three different dot lev¬ 
els all exhibit a sharp current peak initially, because the 
dot to leads coupling is suddenly switched on. For the 
dot level well above the Fermi level, negative value of 
Idisp is found immediately after the initial current peak, 
and indicates that transient charge depletion processes 
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in the center region can also exist in this sudden switch¬ 
ing on case. At the time t fa 27rn/a;o, there are also some 
peak or dip structures at the displacement current, which 
manifest the EPI effect. 


B. Currents driven by an upward step pulse 

Now, we consider another situation with the quantum 
dot initially coupled with left and right leads, and the 
system is in an equilibrium state with zero bias voltage. 
Then a finite bias applied symmetrically between leads 
(/ii = = A^/2) is suddenly turned on at the time 

to = 0, and a transient current is driven through the 
quantum dot. 



FIG. 4. Left, right and total current /L(t), Init) and 7(t) 
after an upward step pulse. The dot level td = 5.0, and 
the electron-phonon coupling constants are: (a) g = 1.0; (b) 
g = 2.0 , respectively. The bias voltage for the pulse g,L = 
—g,R = 2.0. 

In Fig. 4 the time-dependent left, right and total cur¬ 
rent after an upward step pulse are plotted. Both in 
Fig. 4(a) and (b) sharp rising of current amplitude af¬ 
ter turning on the bias voltage are observed. For weak 
EPI case in Fig. 4(a), all the left, right and total cur¬ 
rents have one broad peak initially, then there are some 
smooth current oscillations as the system evolves to the 
steady state. The magnitude of the left current Ihit) 
for electron tunneling into the quantum dot is not equal 
to that of the tunneling out current —lR(t), which indi¬ 
cates that there are charge accumulating and depletion 
processes in the centre region. For the system with strong 
EPI [g = 2), the first peak for the current becomes quite 


sharp, and the current amplitude also decreases drasti¬ 
cally compared with the weak EPI case, which can be 
attributed to the strong renormalization of the tunnel¬ 
ing coupling strength F by EPI. In the time evolution of 
the currents, some periodic structures with a time-period 
of 27r/a;o are clearly observed, and it manifests the EPI 
effect in the system. 

Next we study the total current under different step 
voltage pulses in the strong EPI case. Fig. 5(a) and (b) 
correspond to the systems with different dot levels, re¬ 
spectively. When the dot level id is well above the Fermi 
energy {Er = 0) as shown in Fig. 5(a), signihcant in¬ 
creasing of the steady currents in the long-time limit is 
observed with increasing the bias voltage, whereas for the 
dot level located around the Fermi energy shown in Fig. 
5(b), the steady current values are almost independent 
of the bias voltages. It is easy to understand this behav¬ 
ior since in the later case the dot level is already located 
in the transport window even for a small bias voltage. 
Sharp initial current peaks are observed in both cases of 
dot level locations. Broad current peaks following the 
sharp peaks are exhibited for the system with id = 0, 
indicating a relatively larger current flow in this system 
after switching on the bias voltage. As the bias voltage 
increases, peak or dip structures of the total current aris¬ 
ing from EPI interaction at the time t = 2mT/ujo become 
more significant. 



FIG. 5. The time-dependent total current I{t) for the quan¬ 
tum dot system under different upward step bias voltages. 
The dot level are: (a) id = 5.0; (b) id = 0.0, respectively. 
The EPI constants is taken as g = 2.0. 
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V. CONCLUSION 


We have investigated transient effects in the electron 
tunneling through a SMJ coupled with a vibrational 
mode by using the nonequilibrium GF method and a 
propagation scheme for current matrices. We mainly 
focus on the system with the phonon mode frequency 
Wo ^ r and in the strong EPI region. The total cur¬ 
rent when the tunneling to the leads is suddenly switched 
on gives evident EPI effects, such as the emergence of 
plateau structures with a time period of 27r/wo, which 
reveals the frequency of vibrational mode in the molec¬ 
ular junction. The long time for the system to reach a 
steady state represents a strong renormalization of the 
tunneling coupling between leads and the quantum dot 
by the EPI. As for the system under a upward step bias 


pulse, sharp initial current peaks are observed, but the 
oscillation of the total current is significantly suppressed 
in the strong EPI case. EPI effects are also noticeable 
in the time trace of the total current by changing the 
bias voltage or tuning the energy level of the quantum 
dot. One may expect that the main features of transient 
dynamics of the SMJ shown in this work can be testified 
experimentally using ultrafast spectroscopy techniques. 
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